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^ Abstract. We are interested in the kernel of one-dimensional diffusion equations with con- 

Q tinuous coefficients as evaluated by means of explicit discretization schemes of uniform step 

h > in the limit as h — > 0. We consider both semidiscrete triangulations with continuous 
time and explicit Euler schemes with time step small enough for the method to be stable. We 
H find sharp uniform bounds for the convergence rate as a function of the degree of smoothness 

which we conjecture. The bounds also apply to the time derivative of the kernel and its first 
l l two space derivatives. Our proof is constructive and is based on a new technique of path 

conditioning for Markov chains and a rcnormalization group argument. Convergence rates 
V-^ depend on the degree of smoothness and Holder differentiability of the coefficients. We find 

that the fastest convergence rate is of order 0(h 2 ) and is achieved if the coefficients have 
a bounded second derivative. Otherwise, explicit schemes still converge for any degree of 
4-^ Holder differentiability except that the convergence rate is slower. Holder continuity itself is 

not strictly necessary and can be relaxed by an hypothesis of uniform continuity. 
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^ 1. Introduction and Notations 

& Consider a pair of backward and forward one-dimensional diffusion equations of the form 

(1-1) ~f(x;t) + C°J(x;t)=0, ~g(y; t) = L^gfa t) 

where 

(i.2) C ° = r^ 2 & + ^1- 

and its adjoint formally acts as follows: 

on a test function (p. These equations are defined on the bounded interval A = [— L, L] C K 
where < L < oo. For simplicity, we assume periodic boundary conditions and identify the two 
boundary points ±L with each other. 
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In most of the paper, the coefficients <j(x) and p{x) are assumed to be Holder differentiable. 
More precisely, if a G (0, 1], k € N where N = {0, 1, ...} and the function <j>(x) € C k {A) has k > 1 
continuous derivatives, then one says that <f> is Holder differentiable of order (k, a) if there exists 
a constant c > such that 

(1.4) \^ k \x)-^{y)\<cd{x,y) a 

uniformly for all x,y £ A. In case k = the function is called Holder continuous. The distance 
is defined consistently with the assumed periodic boundary conditions and is given by 

(1.5) d(x, y) = min\x — y — 2Ln\. 

n 

The linear space of Holder continuous or Holder differentiable periodic functions of order (k, a) 
on A is denoted with C k,a (A). We are interested in the case where a 2 € C k ' a (A) and p € C j ^(A) 
with k + a > and j + (3 > 0. 

The hypothesis of Holder continuity can be relaxed slightly by assuming uniform continuity 
instead, i.e. that both u(x) 2 and /j,(x) satisfy a bound of the form 

(1.6) \<f>(x)-<t>(y)\<cp(d(x,y)) 
where p{d) is a non-decreasing function such that lim^o p{d) = 0. 



Let Uq(x, y; t) denote a kernel of equation ( 1.17 1, i.e. a weak solution of the forward equation 



d 

(1.7) —u (x,y;t)=£°*u (x,y;t) 

where the operator £®* acts on the y coordinate and the following initial time condition is 
satisfied: 

(1.8) ]imu (x,y;t) = 8(x- y). 

The kernel uo(x,y;t) formally satisfies also the backward equation 

d 

(1.9) — u (x,y;t) + C°u (x,y;t) = 

where the operator £® acts on the x coordinate. 

We are interested in existence, uniqueness, smoothness and approximation schemes for the 
kernel Uo(x,y;t), its first two space derivatives with respect to the x variable and its first time 
derivative dtUo{x,y\t). As a byproduct of this analysis, we also find conclusions about the 
convergence of £®uq(x, y; t) and L^*uq(x, y; t), as both expressions equal the first time derivative. 

Diffusion equations are one of the single most studied topics in the literature. Existence and 
uniqueness questions for the kernel were address in (Kolmogorov 1931), (Feller 1936), (Hillc 
1948), (Yosida 1951) and (Ito 1957). A classification of all the possible boundary conditions is in 
(Feller 1952). The case of Holder continuous coefficients was resolved in (Philips 1961) based on 
methods in (Friedrichs 1958) and (Lax and Phillips 1960). The hypothesis of Holder continuity 
was relaxed to uniform continuity in (Fabes and Riviere 1966) and (Stroock and Varadhan 1969). 
Strook and Varadhan also introduce a new probabilistic framework where existence is proved by 
reduction to the so-called martingale problem and a compactness argument, thus shifting the 
attention from the kernel itself to the underlying measure space. 

The existence of a weak limit of continuous time Markov chains as h m | was established in 
(Sova 1967) and (Kurtz 1969) by using operator semigroup methods, see also the book (Ethier 
and Kurtz 1986) for a review. Convergence in the semigroup sense takes place if the limit 

(1.10) {Tt<t>){x)= Hm V u hm (x,y;t)4>(y) 

m — >oo, m>n L — ^ 

yeh m ZnA 

exists for all test functions </> € C°°(A), uniformly for all x G A n ,n > 0. A key result is that 
a necessary and sufficient condition for this limit to exist and define a semigroup T t is that 
generators converge also in the same Banach space, i.e. that also the limit 

(1.11) lim £' lm ^ = Cl4> 

hmiO 
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exists in the uniform norm for all test functions cf> € C°°(A). See (Ethier and Kurtz 1986) for a 
precise statement with all the technical conditions and a proof. In (Stroock and Varadhan 1979) 
convergence is reconsidered again by reduction to the martingale problem. 

The problem has also been studied extensively in the numerical analysis literature. Explicit 
and implicit Euler schemes where coefficients are smooth and the data is rough in the sense that 
it belongs to a L 2 space have been considered by several authors. In the case that the Markov 
generator is symmetric and time independent, one can make use of a spectral representation as 
in (Baker et al. 1977) and with greater effort such methods may also be used for more general 
situations, see (Suzuki 1978). In (Luskin and Rannacher 1978), a parabolic duality argument 
is used to show convergence for the standard Galerkin method. (Mingyou and Thomee 1982) 
use a simpler argument based on energy estimates. In (Palencia 1996) one finds convergence 
bounds in maximum norm assuming the initial condition is uniformly bounded and coefficients 
are constant. 

In this article, we revisit this classic theme by considering the problem of obtaining the kernel 
constructively as a limit of increasingly fine triangulations schemes and in assessing the rate of 
convergence with pointwise bounds on the kernel itself. More precisely, let h m — L2~ m ,m S N 
and let A m h m Z n A. Consider the sequence of operators 

(1.12) £- = ^A- + /i(z)V- 

defined on the 2 m+1 -dimensional space of all periodic functions f m : A m — > R, where 

(1.13) V™ f{x) = + 
and 

/-i -i a\ Am ft \ f(x + h m ) + f(x-h m )-2f(x) 

These definitions also apply to the boundary points by periodicity. We assume that m > m 
where mo is the least integer such that 

(1-15) ^ > HX)1 



2hl 2K 



for all m > to and all x G A r< 



Let u m (x, y; t) denote the kernel of equation (1.171, i.e. the solution of the (forward) equation 



(1-16) ^_u m (x,y;t) = C™*u m (x,y;t) 

where the operator C m * acts on the y coordinate and the following initial time condition is 
satisfied: 

(1.17) \imu m (x,y;t) = 5 m (x - y). 

Here, 

. , if x = y mod 2L 

(1-18) 6 m (x-y) = {*™ y 

otherwise. 



Since (1.171 is a finite system of linear ordinary differential equations, the solution exists and is 



unique for all times. The kernel u m (x,y;t) satisfies also the backward equation 

d 

(1.19) 77-u m (x, y; t) + £™u m (x, y,t) = 

at 

where the operator C™ acts on the x coordinate. Using functional calculus notations for the 
exponential of a matrix, we also have that 

(1.20) u m {x,y;t) = exp (tC m )(x, y). 

fl rn 

Our main result can be stated as follows: 
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Theorem 1. Suppose that a 2 6 C k,a and p 6 C j ' and let 
(1.21) 7 = min{2,fc + a, j + (3}. 

Assume that 7 > 0. Then there is a constant c > such that for all m! > m > m Q and all 
x,y G A m the following inequalities hold: 

(i) 

\u m (x,y;t) - u m i(x,y; t)\< cli^ 



(1.22) 
(ii) 



\d t u m {x, y;t) - d t u m -(x 1 y\t)\ 

= \£™*u m (x,y,t) - £™'*u m ,(x,y;t)\ 



l {x,y;t) - C™ u m ,(x,y;t)\< ch? m 



(1.23) 

A version of this theorem under slightly weaker conditions can be formulated as follows: 

Theorem 2. Suppose that a(x) 2 and p(x) are uniformly continuous functions in A. Let the 
function p(d) be non- decreasing and be such that lim^o p(d) = and equation (1.6) holds. Then 
there is a constant c > such that for all m! > m > tuq and all x, y £ A m the following 
inequalities hold: 

(i) 

(1.24) \u m (x,y;t) - u m >{x,y;t)\< cp(h m ) 

(ii) 

\d t u m {x,y;t) - d t u m >(x,y;t)\ 

= \£" y l *u m (x,y;t)-£™'*u m ,(x,y;t)\ 

= \£7 u m{x,y;t) - C™ u m >(x,y;t)\< cp(h m ). 

(1.25) 

Next, we consider the case where also time is discretized and prove the following result: 



Theorem 3. Suppose that a 2 and p satisfy equations of the form (1.6) with a non- decreasing 
function p(d) such that lim^o p{d) = 0- Consider the discretized kernel 



(1.26) 



u s r l(x,y;t)=h- 1 (l + StjC m )^(x,y;t). 



where C m is the operator in (1.12) and St m is so small that 
(1.27) min 1 + St m C m (x, x) > 

x£A m 

Assume that boundary conditions are periodic and that the ratio ^ — N is an integer. Then 
here is a constant c > such that the following bounds hold for all m > mo and all x, y € A m : 

(i) 

\u m (x,y;t) - u%(x,y;t)\< ch 2 m 



(1.28) 
(ii) 



u s ^{x,y;t + St)-u s ^x,y;t) 
St 



d t u m (x,y;t) 

i(x,y,t) - £ 
\£Tu m (x, y, t) - CTut(x, y;t)\< ch 2 



= \£™*u m (x,yt) - £™*u%(x,y;t)\ 



(1.29) 
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The paper is organized as follows. In Section 2 we consider the case of Brownian motion and 
review a result in (Albanese and Mijatovic 2006) which establishes the theorems above in this 
simple particular case where Fourier analysis in the space direction can be used to carry out a 
precise calculation. In Section 3, we consider the case of a diffusion where both the volatility 
and the drift have two bounded derivatives. In this case, we make use of time-homogeneity and 
carry out a Fourier transform in the time direction after path conditioning. In Section 4, we 
extend the derivation to the case of non-smooth coefficients. Section 5 is dedicated to the case 
where time is discretized and we prove Theorem [3] 



2. Constant Coefficients 

In this Section, we prove Theorem 1 in the special case where the volatility and the drift 
coefficients are constant, i.e. 

1 



(2.1) 



I- x = /iv x + -a A x . 



It suffices to consider the case ml = 



1. Let B rn be the Brillouin zone defined as follows: 



(2.2) 

Let T m : i 2 {A Tl 
(2.3) 



B m — I — 



kir 



,fc = 0, ..2 m -l 



£ 2 (B m ) be the Fourier transform operator defined so that: 



f{p)=T m {f){ P ) = K 



me 



for all p £ B m . The inverse Fourier transform is given by 



(2.4) 



^•- 1 (/)(x) = ^ ; hp) 



p£B„ 



The Fourier transformed generator is diagonal and is given by the operator of multiplication 



by 
(2.5) 

We have 
(2.6) 



sin h m p 



> cos h m p - 1 
hi • 



u m (x,y;t) 



2L ^ 

p£B 



\^ e t^ m (p) e ip{y-x) 



Using this Fourier series representation, we find 
(^} y ;t) - u m+1 (x,y;t)\ 



< 



1 

2L 



E 

p£B m 



j e ip(y-x) 


1 




+ 2L 



E 

p£B m+1 \B, 



e ti m+1 (p) e ip(y-x) 



(2.7) 

Let 

(2.8) 



K m = 



| log h 



m+l | 



If h m is small enough, i.e. if toq is sufficiently large, we have that 



1 

2L 
(2.9) 



E 

P eB m .\p\>K rr 



e tt™{p) e ip(y-x) 



< 



- E 

pEB m ,\p\>K„ 



e mi™( P )) < cexp ( ta 



j2 cos h m K 7n 1 

h 2 
1 1 ... 



< chi 
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where 3?(a) denotes the real part of a G C and c denotes a generic constant. Similarly 



1 

2L 



E 

p£B m+1 ,\p\>K„ 



e a m+1 (p) e w(y-x) 



< 



1 

2L 



E 



< cexp to 



peB m ,\ P \>K m 
cos h m+ \K — 1 



iSR(^ +1 (p)) 



(2.10) 
Since 

(2.11) 

and 

(2.12) 



2h ~ 2 P 



(2ft)s 



we find that if \p\< ^ then 
(2.13) 

Moreover, since 
(2.14) 



\r(p)-r + \p)\<^h 2 \p\ 3 +^hY. 



1 2 cos hp — 1 1 



ft 



we conclude that in case |p|< h 1 \J\, the following inequality holds: 



(2.15) 

Hence, if m is large enough, we find 



cos hp — 1 In 

< — v 

h ~ V 



1 

2L 



E 

peB m ,\ P \<K 



,u m (p) _ (p) ) Jp{v-x) 



< 



(2.16) 



< 



L E 

2L ^ 

pS-B m ,|p 

1 E 



pe-B™,|p|<^ 



e 4f 



^IpI^+t^^p 4 _ i 



4 <br+^> 4 ) < chi 



p£B m ,\p\<K 

for some constant c > independent of m. This concludes the proof of convergence for 
kernel in the special case of constant coefficients. 
To estimate the first derivative, notice that 



(2.17) 
and 



(z, 2/;*) 



< 



+ 



i 

2L 
1 

2i 



E 

pes™ 



; « m (p) 



sinp/i 



m _ e te m+1 ( P ) 



suiph 
h m +i 



±1 | e ip(y-x) 



E 

pG-B m+ i\B, 



e ti m+1 (p) e ip(y-x) 



(2.18) 
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Let 
(2.19) 

If h m is small enough, we have that 
1 

2L 



K — 2 



|logfem+ll 



E 

p£B m ,\p\>K„ 



sin ph„ 

h,: 



e U m (p) °"^»m c ip(y-x) 



< 



2L 



E 



p£B m ,\ P \>I<n 



hr, 



< c 



sin _?£T/i, 



exp ta 



, cos Xft- In — 1 

' 'mi 



< chi 



(2.20) 

where c denotes a generic constant. Similarly 

1 

2L 



(2.21) 



E 

pGB m+ i,b|>-RT„ 



e « m+1 (p) e v(a-K) 



< ch J 



If 77i is large enough, we also find 



1 

2L 



E 

peB m ,\p\<K r , 

<- E 



sin p/7 



m e u m (p) 



sin pft,,, 



h m +i 



+1 „« m+1 (p) \ p ip{y-x) 



sin ph 7n 









peB m ,| P |</Y 

sinp/i„ l+ i sinp/i 



h m +\ 



< chi 



(2.22) 



for some constant c > independent of m. This concludes the proof of the bound of the first 
derivative. The second derivative can be derived in a similar way. 

Finally, consider the following Fourier representation for the discretized kernel 



.St 



(x,y;t) = j ]T (l + 5tl™{p)Y~ e 



peB„ 



ip(y-x) 



(2.23) 

Consider the formula 
(2.24) 

and let's represent the difference between the discrete and continuous time kernels as follows: 
\u m (x,y;t) - u & ^(x,y\t)\ 



1 + 8tt m {p) ) = exp ( t log (l + £™(p)) 



< 



1 

2L 



peB„ 



]T ( exp (ti m (p)) - exp ( y log (1 + 8t( m {p)) ) e^"-) 



< 



2L 



1 y 

eB m ,p<K m 

1 E 



e 4 ' 



St 



exp ( - log (1 + 5U m {p)) - U m (p) ) - 1 



p£B m ,p>K„ 



exp (ti m (p)) 



- E 



p£B m ,p>K„ 



exp(-log(l + <Sir>)) 



(2.25) 



where K m is chosen as in (2.8). The very same bounds above lead to the conclusion that this 
difference is < ch"L . 
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Re(o) 



Figure 1. Contour of integration for the integral in (3.521. C+ is the countour 
joining the point D to the points E, A, B. C_ is the countour joining the point 
B to C to D. 



3. Smooth Coefficients 

In this section, we prove Theorem 1 in the case where the drift and volatility are both of class 
C^ 3,0 ), i.e. they depend smoothly on the space coordinate but not on the time coordinate. 
Let us introduce the following two constants characterizing the volatility function: 



(3.1) S = inf a(x), Si = sup \/cr(x) 2 + h m \fi(x)\. 

*£A m xGA m 

and let 

(3.2) M = sup |ju(x)|. 

xeA m 

Since our interval is bounded, we have that So > and Si, M < oo. 

A symbolic path 7 = {70, 71, 72, ■•■•} is an infinite sequence of sites in A m such that "fj ^= Jj-i 
for all j = 1, .... Let T m be the set of all symbolic paths in A m . The kernel of the diffusion 
process admits the following representation in terms of a summation over symbolic paths 

^ 00 

(3.3) u m (x,y;t) =—J2^ q E W m (j,q,t) 

q ~ 7 S T m : 70 = x,j q = y 
l7j-7j-i|= 1 V j > 1 

where 
(3.4) 



f dsx f ds 2 ... f ds q e^- s ^ c ^'^ ]J (e^-^^^hCmi^nj+i)) 

flm Jo Jsi Js q -i 1=n V / 



with sq = 0. 
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Let us introduce the following Green's function: 



(3.5) 



G m (x,y;u) 



u m (x,y; t)e~ 



1 



The propagator can be expressed as the following contour integral 



— (x,y). 



(3.6) 



u m (x,y;t) 



2^ 



G m (x,y]U))e l 



dui 
2n 



G m (x,y;cj)e 1 ' 



Here, C+ is the contour joining the point D to the points E, A, B in Fig. [T] while C_ is the 
contour joining the point B to C to D. By design, each point oj on the upper path C+ is 
separated from the spectrum of C. 

Lemma 1. For m sufficiently large, there is a constant c > such that 



(3.7) 



2vr 



G m (x,y;u))e 1, 



< cbf. 



Proof. The proof is based on the geometric series expansion 

1 . ,^ 1 



(3.8) 



G m (u>) = h Ti 



= h: 



whose convergence for uj € C+ can be established by means of a Kato-Rellich relative bound, see 
(Kato 1966). More precisely, for any a > 0, one can find a (3 > such that the operators V m 
and A m satisfy the following relative bound estimate: 



(3.9) 



||V m /|| 2 <a||A m /|| 2 +/3||/|| ; 



for all periodic functions / and all m > mo. This bound can be derived by observing that V m 
and A m can be diagonalized simultaneously by a Fourier transform, as done in the previous 
section, and by observing that for any a > 0, one can find a [3 > such that 



(3.10) 



sin h m p 



< a 



cos h m p — 1 



h 2 



for all m > mo and all p € B m . 

Under the same conditions, we also have that 



(3.11) 
Hence 
(3.12) 



< 



2Ma 



Y 2 



1 



-a 2 A m f 



icr 2 A T ' 



/ 



< 



2Ma 
y-2 



-a 2 A r 



\o 2 * 



■0 



f 



no 



■0 



\o 2 * 



-f 



< 1 



where the last inequality holds if lj G C+, if a is chosen sufficiently small and if m is large enough. 
In this case, the geometric series expansion converges in (3.8 1 converges in L? operator norm. 
The uniform norm of the kernel \G m (x,y;uj)\ is pointwise bounded from above by h^. 



Since the points B and D have imaginary part equal at height 4 



the integral over the 



contour C+ converges also and is bounded from above by ch m in uniform norm. 



□ 



Lemma 2. If q > 

(3.13) 



Eft 



we have that 



W m { 1: q-t) < 
Proof. Let us define the function 
(3.14) 



exp 



m = 2fk 6 1{t -° ] 
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where l(i > 0) is the characteristic function of R+. We have that 

(3.15) w m (7,<z;*) <4>* q (t) 

where <j)* q is the q— th convolution power, i.e. the q— fold convolution product of the function 
by itself. The Fourier transform of </>(i) is given by 

y^2 i-oo Egt V 2 

(3.16) = e-^-^dt- El 



2KL ./n " 2iojh m + Eg ' 



The convolution power is given by the following inverse Fourier transform: 



(3,7) ^ W .Jf *MV"'=^j U' + ^fj *'"V 

Introducing the new variable z = 1 + 2 'y2 m , the integral can be recast as follows 

y2-2<7y2<; / y2y- \ 

where C# is the contour in Fig. [2] Using the residue theorem and noticing that the only pole of 
the integrand is at z = 0, we find 

1 /£?t\ 9 /-Eg* 



<3 ' 19) ^> = (^ U4J ^ V 2^ 

Making use of Stirling's formula q\ ss A/^g'+s e~ 9 , we find 

(3.20, » W^exp (-f| + s .„ g f| + ,(1 - log,, 



2 f 

If logq > log 2775- + 2, then we arrive at the bound in (3.13) 



□ 

It suffices to consider the case m! — m+ 1 for all values of m above a fixed threshold. In fact, 
given this particular case, the general statement can be derived with an iterative argument. To 
this end, we introduce a renormalization group transformation based on the notion of decorating 
path. 

Definition 1. (Decorating Paths.) Let m > m and let 7 = {yo,yi,U2, ■■■■} be a symbolic 
sequence in T m . A decorating path around 7 is defined as a symbolic sequence 7' = {yo, y[,y'2, ■■■■} 
with y[ € h m+ iL containing the sequence 7 as a subset and such that if y'j = yi and y' k = g/j+i, 
then all elements y' n with j < n < k are such that \y' n — y'j\< h m+ \. Let T> m+ i{^) be the set of 
all decorating sequences around 7. The decorated weights are defined as follows: 

00 

(3.21) W m ( 7 ,q;t) = Y^ E W m+1 (i ,q' ;t). 

«-« 7 '6P m+ i(7) 

Finally, let us introduce also the following Fourier transform: 

roc A /"OO 

(3.22) W m {j,q;uj)= VF m ( 7 , q; 0, t)e iut dt, W m (j,q;oj)= W m (j,q;t)e iuJt dt. 

Jo Jo 

Definition 2. (Notations.) Ln the following, we set h — h m+ i so that h m — 2h. We also 
use the Landau notation 0{h n ) to indicate a function f(h) such that h~ n f(h) is bounded in a 
neighborhood o/(0). 

Lemma 3. Let A m and let C_ be an integration contour as in Fig. [JJ Then 

, dto 



(3.23) 
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Figure 2. Contour of integration Cr for the integral in (3.18) 



(2W m { 7 ,q;uj)- W m (7, 



Proof. We have that 
(3.24) 

2G m+1 {x,y;u) - G m (x,y;u) = ^E 2 ~ 9 X! 

9=1 7 € r m : 7o = x,J q = V 
l7i-7j-i|= IVj > 1 

The number of paths over which the summation is extended is 

(3.25) N(j,q;x,y) = (({7 € T m : 7o = 21,7, = y, \jj - Tj-i!= IVj > 1} = f « + fc) 

where fc = Applying Stirling's formula we find 



(3.26) 

Hence 



^7 < 2« 



2G m+ i(x, y; w) - G m (x, y; w) e 



OO 

C x > /l 



max 

g=l ■ 9 7 G T m : 70 = x, 7, = y 
l7i-7j-i|= IVj > 1 



, duj 
2tt 



2W m ( 7 ,q;u;)-W m ( 1 ,q;u ) )je iut ^- 



(3.27) 
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for some constant c w W ^ > 0. It suffices to extend the summation over q only up to 



(3.28) 



e 2 £ 2 i 



ymax- ^ ■ 

To resum beyond this threshold, one can use the previous lemma. More precisely, we have that 



1( 



2G m+1 (x,y;uj) - G m {x,y;uj) e 



Aujt 



dw 
2n 



< 



h 



2^ 



max 

g,7 € T m : 7 = x,j q = y 
l7i-7j-i|=lVi>l 

(3.29) 

Let v(x) = <r(x) 2 . To evaluate the resummed weight function, let us form the matrix 





(3.30) C{x;h) = 

and decompose it as follows: 



v{x+h) 

TP 

lb? 2h 




v(x+h) ix{x+h) 
~^SP 7~2K 

v(x) 

~~TP~ 
v(x—h) _|_ /j,(x—h) 
~^TP 1 2h 



v(x) fJ-(x) 
~2~TP~ W 
v(x—h) 

TP 



(3.31) 
where 

(3.32) 



£(x; h) = p£ (x) + ^i(x) + £ 2 (x) + h£ 3 (x) + 0(h 2 ). 



-v(x) \ v {x) 
C {x) = I \v{x) -v(x) \v(x) I , 



\v{x) -v(x) 



(3.33) 

(3.34) 

and 

(3.35) 



-v\x) 

d{x) = I y(x) 



\v'(x) - \n(x) 


2 \ X ) + 2 ' 





-V(a;) + hl*(x) v'{x) 



C 2 {x) = 



■\v"{x) 





-\v'"{x) 




\v"{x) y{x) 







i„,« 

4 







!„,»/' 



v"(x)-^'(x) -\v"{x) 









-i«'"(x) + i M "(x) § V '"(x) 








Let us introduce the sign variable r = ±1, the functions 

(3.36) Mt,x,r) = 2£ m (x,x + 2Th)e tCm< - x ' x h(t > 0) 

(3.37) <l>i{t,x,T) = 2£ m+1 (x + Th 1 x + 2Th)e t£(x ' h) (x 1 x + Th)l(t > 0) 
and their Fourier transforms 



<j>o{u),x,T) 
(3.38) 



v(x) n{x) \ ( v(x) 



2h / V Ah 2 



v(x) , fi(x) + v'{x) , v"{x)+fj,'{x) , ^'"(x) , ^"(x) 
1? 



+ r 



+ 



+ 



6 + 2 



r/ l + 0(/i 2 ) 



< x| (— £(x; h) + iuj) 1 \x + rh > . 
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where 

M / s rA 

(3.39) \x>= 1 , and \x + rh >= 

W \«r,-l 

We also require the functions 

(3.40) ijj (t,x) = e tCml - x ' x h(t>Q), ip 1 (t,x) = e t£{v:h) {x,x)l(t>0) 
and the corresponding Fourier transforms 

(3.41) fo(w,x)=(?& + iw\ , Miu,x)=<x\(-£(x;h) + iuy 1 \x> . 
If 7 is a symbolic sequence, then 

9-1 ^ 

(3-42) W m (-y,q;uj) = ip (ui,j q ) J| <Mw; 7,, sgn(j j+ i - jj)) 

3=0 
9-1 ^ 

(3-43) # m (7,g;w) = ^(w, 7,) <£i(w; 7,-, sgn( 7j+ i - 7,-)). 

j=0 

Let us estimate the difference between the functions </>i(w, x, r) and ^(w, x, r) assuming that 
w is in the contour C_ in Fig. [2j Retaining only terms up to order up to 0(/i 3 ), we find 

(3.44) 

A lengthy but straightforward calculation which is best carried out using a symbolic manipulation 
program, gives 

1 1 \ 2a(x)Th Aiwti 2 r .iuirh 3 

M», x,r) = l+ - ^ - - - (*)] ^ 



1 a, t 2u4 

+ r{x) ■ h 3 r + iuh A p{x) — - + 0(h 5 ) 

v\xy 



(3.45) 
where 



r ( x ) = r, 1 ^ \t*"( x ) v ( x ) ~ 4w(x) 3 + 2v'(x)fi(x) 2 — 2v'(x)v(x)fi'(x) 
Zv(x)- 3 

(2fj,(x)fj/(x) + v"{x)v{x) - 2v'(x) 2 ) f i(x)] . 

p(x) = — ^ \4:fi(x) 2 - 2v'(x)u(x) + 4v(x)u'(x) + v"(x)v(x) - 2v'(x) 2 } . 
v{x) i L ' J 

(3.46) 

We have that 
9-1 / 

^2 ( lo g 0o (w; 7j,sgn(7j + i -7j))-log<£i(w;7j,sgn(7 J - + i-7 i )) 

= £ + K7,))^sgn( 7j+ i - 7,) + (^||b||oo+2|^| 2 || V - 2 |U)0(/ i 4 <z) 

= z^ 2 log f^j) + h 2 {R( lq ) - i?( 7o )) + (M||p|U+2M 2 ||OL)0(^) 

(3.47) 



14 



CLAUDIO ALBANESE 



where R(x) is a primitive of r{x), i.e. 
(3.48) R(x) 



r{z)dz. 



We conclude that there is a constant c > such that 

,q-\ q-l 

(3.49) 

V 'j=o J=0 



/ ( TT <fo(w;7j,sgn(7j+i - 7i)) - TT ^i(w;7j,sgn(7j+i - 7j)) 



2tt 



< c/i 2 . 



for all q < q max . Here we use the decay of e lult in the upper half of the complex u> plane to offset 
the u> dependencies in the integrand. Similar calculations lead to the following expansions: 



4h 2 



2h 2 



1 



(3.50) ^ O (u,y) = ^+O(uh i ) ) ^( Wl y)=— + OK) = ^,!/) + OK)- 
v(y) v(y) 2 



Since q < ch 2 and u < \\ogh\, we find 



(3.51) 



c_ 



2G m+ i(x, y; u) - G m (x, y; u>) e 



, dui 
2k 



< c 



: -h A < ch 2 . 



This completes the proof of the Lemma and of the Theorem. 

By differentiating with respect to time in equation |3.52[ we find that 



□ 



(3.52) 



d_ 
Of 



u m {x,y;t) 



du 
2^ 



itoG m (x,y;uj)e 1 ' 



du> 
2^ 



iujG m (x,y;uj)e' 1 ' 



All the derivations above carry through and we conclude that 



(3.53) 

and also 
(3.54) 



iujG m (x,y;u)e l 



2tt 



iu 2G m+ i(x, y; u) - G m (x, y;u) e 



< ch* 



, du> 
2?r 



< c^"^h A < ch 2 . 
h 



Hence the first time derivatives of the kernel satisfy the same Cauchy convergence condition as 
the kernel itself. 

4. Lesser Smooth Coefficients 

In this section we assume coefficients are either Holder continuous or obey the conditions in 
Theorem 2. 

Lemma 4. Let f{x) be a continuous junction in [-L, L] satisfying periodic boundary conditions. 
Then, for all h > 0, we have that 



(4.1) 
(4.2) 



f{x + h m ) = f{x) + h m V™f(x) + — A™f(x) 

f{x - hm) = f{x) - h m V m fix) + yA™/(x). 



This is the result of a simple calculation, which is however useful as it allows one to extend 
the derivation in the previous section by making the following replacements: 



(4.3) 
(4.4) 

In fact, 
(4.5) 



v'ix) V>(z), v"ix) -> A m vix), v"'{x) -> 
l/(x) - V>(z), n"ix) -> A™nix). 



C(x; h) = -p^Coix) + \Ciix) + C 2 ix) + hC^x) 
h z h 
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without any 0(h 3 ) corrections as long as one re-defines the matrices on the right hand side as 
follows: 



(4.6) 



'— v(x) \v{x) 







£o(x) = \v{x) -v{x) \v{x) , 



\v(x) 



-v(x) / 



(4.7) 



Ci(x) 



-V?v(x) \V™v{x) - \ii{x) 



\n{x) 






-\n{x) 
\V™v{x) + y(x) V?v(x), 



(4.8) 

and 

(4.9) 



C 2 (x) = 



■\A™v{x) 



Cz{x) 



\*gv{x) - §V™/x(*) 


\A™v(x) - iV?M*) 






\^v(x) ; 



<o 4A>(x) o^ 





All derivations in the previous section go through formally unchanged and one arrives at the 
following expressions 



0o(w, x,t) = 1+ ~-~Y\ -— - 8fi(x) T -r r + 0(h 5 ). 



v(x) v(x) 



v(x) 2 w(x) 2 



(4.10) 
and 



x, t) = 1 + ^ - -j^ - [8»(x) - VXz)] 



u(x) w(x) 



+ r{x) ■ h 3 r + iuuh 4 p(x) - '"7 " + 0(/i 5 ) 



Uco 2 h 4 
v{x) 2 



(4.11) 
where 
r(x) 



p(x) 
(4.12) 



^3 [AJXzM*) - V(x) 3 + 2V;^(x) A1 (x) 2 - 2V>(2>(x)VX*) 

- (2 M (x)V" M (x) + AX,)^) - 2V>(x) 2 ) M (x)] . 
■ [Mx) 2 - 2V™v(xMx) + 4v(x)V?»(x) + A™v(x)v(x) - 2V>(x) 2 ] . 



v(x) : 



We have that 
9-1 



( log0o(w;7 J ,sgn( 7j+ i - 7j -)) - log <£i(u;; 7,-, sgn(7 j+ i - lo )) 



j=0 
< 



" 1 ^V%tt( 7j -) 



E 

3=0 



+ (|c|||p|| co +2|c| 2 || U - 2 || co )0(/ l 4 g). 



< 2M 2 sup 



iwV!?t>(a;) 
^—t^ +r{x) 



v(x) 2 



+ (^||H| co +2|c| 2 ||«- 2 || co )0(/ i 4 g ) < ch\ 



(4.13) 
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where in the last step we made use of the Holder continuity assumptions of Theorem 1. The 
other bounds staying the same, we arrive at 



(4.14) 



2G m+1 (x, y; u) - G m (x, y; w) ) e wt <U ' 



2tt 

Under the weaker assumption of Theorem 2, the bound that applies is instead 



(4.15) 



2G m+ i(x,y;uj) - G m (x,y;uj) )e' ( 



< c^--h 2 p(h) < cp(h). 



Similar bounds also extend to the case of the first time derivative, since multiplication by a 
factor iu> inside of the contour integral is immaterial as far as establishing a bound of this sort 
is concerned. This completes the proof of Theorem 1 and Theorem 2. 



5. Explicit Euler Scheme 

In this section we prove Theorem [3] A Dyson expansion can also be obtained for the time- 
discretized kernel and has the form 



oo N N N 

«S(».»it)TE E EE- E 

9=1 7er m :7o=£C,7<j=ly fci=l fe 2 =fci+l k q =k q -i+l 

/ \ fei— 1 g / \ fej+i— kj— 1 

(5.1) ^l + ^£ m ( 7o ,7o)J («5t) 9 J] J C m (7,_ 1 ,7,)^l + ^£ m (7,, 7 ,)J 

where t q +i — t and k q+ \ — N. In this case, the propagator can be expressed through a Fourier 
integral as follows: 

(5-2) u%(yi,y 2 ;t) = /" G^y^e^ 



2tt 



where 



(5.3) G%(yi,y 2 ;to) - StJ^u^yxJS^e-^. 

3=0 

The propagator can also be represented as the limit 



(5-4) u m (yi,y 2 ;t) = Km / G m (yi,y 2 ;uj)e 

H^oc J Cr ITT 

where Ch is the contour in Fig. [3] This is due to the fact that the integral along the segments BC 
and DA are the negative of each other, while the integral over CD tends to zero exponentially 
fast as 9(w) — > oo, where 3(w) is the imaginary part of u. Using Cauchy's theorem, the contour 
in Fig. [3] can be deformed into the contour in Fig. [T] To estimate the discrepancy between the 
time-discretized kernel and the continuous time one, one can thus compare the Green's function 
along such contour. Again, the only arc that requires detailed attention is the arc BCD, as the 
integral over rest of the contour of integration can be bounded from above as in the previous 
section. 

Let h = h m and let us introduce the two functions 

(5.5) <M*,£,t) = 2C rn (x,x + Th)e tCm(x ' x) l{t > 0), 

(5.6) 4> St {j,x,T) = 2C m (x,x + rh)(l + 8tC m {x,x)) 3 . 
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D 



-ir/St 



lm(ra) 



c/h 



iH+w/6t 
C 



Tf/St 



Re(co) 



Figure 3. Contour of integration for the integral in (5.4) 



and the corresponding Fourier transforms 

, > ~ , . Z" 00 , , , i,,tduj ( v(x) u(x)\ ( v(x) 

(5.7) 0o(w,x,r)= / <^o(*,a;,T)e- ,u *— = I 44 + t4-^ 1 I 44 ■+ 



2tt 



/i 2 



/; 2 



(5.8) M", *, r) = g *, r)e-*-i« = + r ^) (> 5t - 1 + St^ " 



We have that 



(5.9) 



.' v(x) K x )\ ( ■ v ( x ) 



0o(w, r) + T^^St + 0{h 2 5t 2 ). = Muj, x, T ) + 0(h 4 ), 
2v(x) 



where the last step uses the fact that St = 0(h 2 ). 
Let us also introduce the functions 



(5.10) Mt, x ,r) = e tc ^ x h(t>0), ^ 5t (j,x : T) = J2 (l + 6tC m (x,x)) 



and the corresponding Fourier transforms 

//Mr! \ 

(5.11) 1po(w,X,T) 



v(x) 



h 2 



fc=l 



1> 6t (w,x,T)= ( e iuSt -l + St 



v(x) 



Again we find that 
(5.12) 



4>q{u,x,t) = 4> St (LJ,x,r) + 0(h 4 ). 
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If 7 is a symbolic sequence, then let us set 

9-1 ^ 

(5-13) W m (j,q;uj) = -0 o (w,7 ? ) J| (f>o(u; Jj, sgn(7 i+ i -7,*)) 

3=0 

9-1 ^ 

(5-14) W%(-y,q;uj) = 7,) Y[ 7/ ^gn^+i - 7j))- 

i=o 

We have that 
(5.15) 

G^(x,y; W )-G m (x, 2 ;; W )=-^2^ ]T (w£*( 7> g; w) - W m (-y, q; w)) . 

9=1 7 G r m : 7o = x,Jg = y 
\lj -7i-i|= IVj > 1 

The integration over the contour in Fig. [T] can again be split into an integration over the 
countour C_ and an integration over C+. The integral over C+ can be bounded from above 
thanks to Lemma [T] Furthermore, we have that 

< ch~ 1 ^/q max max 

q,je T m : 70 = x, j q = y 

hj -7i-i|= IVj > 1 

(5.16) < ch 2 

To bound the time derivative, we have to consider 

(5.17) 

But, since St = 0(h 2 ), also this difference is 0(h 2 ). 

6. Conclusions 

We obtained bounds on convergence rates for explicit discretization schemes to the kernel of 
one-dimensional diffusion equations with continuous coefficients. We consider both semidiscrete 
triangulations with continuous time and explicit Euler schemes with time step small enough for 
the method to be stable. The proof is constructive and based on a new technique of path condi- 
tioning for Markov chains and a renormalization group argument. Convergence rates depend on 
the degree of smoothness and Holder differentiability of the coefficients. The method is of more 
general applicability and will be extended in future work. 
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